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The sensitivity of inflationary spectra to initial conditions is addressed in the context of a phe- 
nomenological model that breaks Lorentz invariance by dissipative effects above some threshold 
energy A. These effects are obtained dynamically by coupling the fluctuation modes to extra de- 
grees of freedom which are unobservable below A. Because of the strong dissipative effects in the 
early propagation, only the state of the extra degrees of freedom is relevant for the power spectrum. 
If this state is the ground state, and if A is much larger than the Hubble scale H, the standard 
spectrum is recovered. Using analytical and numerical methods, we calculate the modifications for 
a large class of dissipative models. For all of these, we show that the leading modification (in an 
expansion in H/A) is linear in the decay rate evaluated at horizon exit, and that high frequency 
superimposed oscillations are not generated. The modification is negative when the decay rate 
decreases slower than the cube of H, which means that there is a loss of power on the largest scales. 



I. INTRODUCTION 

Today's picture of cosmological evolution assumes that 
all large scale structures we observe today developed from 
primordial fluctuations on top of a homogeneous and 
isotropic state of the early universe. These primordial 
fluctuations naturally arise in the context of inflation. As 
a consequence of the accelerated expansion, short wave- 
length vacuum fluctuations are amplified as they exit the 
horizon scale. After the end of inflation, these fluctua- 
tions re-enter the horizon and eventually undergo grav- 
itational collapse. Depending on the total number of e- 
folds of inflation, the structures we observe today may 
originate from fluctuations with extremely small initial 
wavelengths as defined in the homogeneous frame. In 
fact, unless we fine-tune the number of e-folds, the rele- 
vant scales were all well beyond the Planck scale at the 
onset of inflation [l[ . 

Inflation, therefore, effectively acts as a spacetime mi- 
croscope, offering the tempting opportunity to probe very 
high energies by lookingfor signatures in the primordial 
perturbation spectrum 2, 3]. Let us denote the scale at 
which new physics becomes important by A. Depend- 
ing on the nature of the dominant new physics (which 
may or may not be of gravitational origin), A might be 
the Planck scale, the string scale, or below. Demanding 
that the theory yields the usual results in the infrared 
generically gives rise to a suppression of the corrections 
by some power of H p /A, where H p is the Hubble scale 
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at the time when the mode p under consideration left 
the horizon. There are thus two ways of detecting the 
new physics in the perturbation spectrum. First, if the 
correction contains a sharply defined phase factor which 
is a function of H p / A, an oscillatory feature extending 
over a wide range of the power spectrum may be pro- 
duced 0, H, @. Second, even if oscillations are absent, 
there is still a possibility that the largest scales we ob- 
serve today correspond to sufficiently large H p /A that 
the new physics lead to a distinctive suppression or en- 
hancement of the low-? power spectrum at a detectable 
level. While we do not find effects of the first kind in our 
model, we conjecture that it generically predicts a large 
scale modification of the power spectrum. 

In the absence of clear predictions from a fundamen- 
tal theory, several types of phenomenological approaches 
have been proposed. If Lorentz invariance is retained at 
all energies, see e.g. Q, there is no possibility to introduce 
deviations from the standard relativistic propagation and 
thus very little hope to get any signatures. It is therefore 
interesting to consider the breaking of Lorentz invariance 
in the ultraviolet (UV) sector, and in fact this was done 
in essentially all approaches. 

In the simplest of these models, scalar or tensor per- 
turbation modes are prescribed to be created in their 
adiabatic vacuum at some fixed initial time, where the 
initial conditions are specified either on a spacelike sur- 
face common to all modes or when the physical mo- 
mentum P = p/a of each mode with comoving momen- 
tum p satisfies P = A 0, d, @|- These models predict 
a spectrum with superimposed oscillations whose ampli- 
tude is a power of H p /A which depends on the degree 
of non-adiabaticity of the initial state. The validity of 
this conclusion was questioned in 0] on the basis that 
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the modulation of the corrections artificially follows from 
the sharp, and thus non-adiabatic, character of imposing 
the initial state at a given instant. 

In another approach [2j, |3[ , Lorentz invariance is broken 
by introducing deviations from the relativistic dispersion 
relation above a UV scale A: 




n+l N 



(1) 



where f2 and P are the proper frequency and the proper 
momentum as measured in the preferred frame which is 
assumed to coincide with the cosmological frame [24| ■ It 
has been understood that the standard predictions are 
robust [12j, i.e. the modifications of the spectra scale as 
a power of H p /A, provided the initial state is the asymp- 
totic vacuum and the modes evolve adiabatically. Under 
these conditions, dispersive models generically predict no 
superimposed oscillations [HI]. 

The alternative possibility that Lorentz invariance is 
broken by dissipative effects has received much less at- 
tention so far. To be realized while maintaining unitarity, 
one must introduce additional, unobservable degrees of 
freedom, hereafter called which couple to the observ- 
able field 4> m the UV sector In this paper we aim 
to compute the modifications of the spectrum induced by 
such dissipative effects. 

In order to obtain a local equation for the effective 
propagation of <f> after tracing out a simple class of 
models for the propagation of 4" and its interaction with 
<j) is analyzed in detail. This class is characterized by the 
"decay rate" T of the <j) modes in the preferred frame, 
which - in analogy to eq. {!]) - we parameterize by 
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If H/A -C 1, we argue in Sec. Ill Fl that the power spec- 
trum of a wide range of dissipative models can be effec- 
tively described by a simplified model characterized by 
such a decay rate. 

The paper is organized as follows. The model is 
presented in Sec. Hi] After introducing the settings in 
Sec. Ill Al we derive the effective equation of motion of 
the 4> mode in Sec. Ill Bl In Sec. Ill Cl the power spectrum 
is presented in terms of a double integral of a noise ker- 
nel governed by \P and the retarded Green function of 4>. 
From an approximate expression of the Green function 
fSec. HII A~]l . we derive analytic expressions for the power 
spectrum in Sec. IIIIBI and IIII CI at zero and high tem- 
peratures, respectively. The numerical scheme and the 
results are presented in Sec. IIVI 



II. MODEL AND POWER SPECTRA 
A. The model 

In slow-roll inflation, the background is a flat Fried- 
mann universe with the usual Friedmann-Lemaitrc- 
Robertson- Walker metric 



ds 2 = -dt 2 



(t) dx 2 



WK+dx 2 ) , (3) 



and the variation of H = <9 t In a is governed by the slow- 
roll parameter e — —dtH/H 2 <C 1. Since the back- 
ground is homogeneous, the fields decompose into Fourier 
modes labeled by the comoving wave vector p. More- 
over, when the cosmological and preferred frames coin- 
cide (see [ll|), the action splits into disconnected sectors, 
S = fd 3 P S(p). 

We consider a scalar field cj). The action for the rescaled 
mode $r) = aS-o is 



where the conformal frequency is given by 



(v) = p 2 



f 



(4) 



(5) 



For the scalar field <fi, one has f — a. The same is true 
for tensor modes, whereas for density perturbation modes 
/ = i/ea. Given the simplicity of these substitutions, we 
limit our discussion to the scalar field in this work. 

The power spectrum of <f> is given by the Fourier trans- 
form of the two-point correlation function at equal times 



P P (l) 



d 3 x 



ipx 



(0 (ry, X) 0(77,0)) 



(2tt) 3 

^ jdVUM^VW}), (6) 



2a 2 (77) 



where ( • ) and { • , • } denote the quantum expectation 
value and the anticommutator, respectively. 

In the Bunch-Davies vacuum (the adiabatic vacuum 
for P = p/a — > 00), the power spectrum is simply 



1 



a 2 (77) 



(7) 



Here, $p denotes the unit Wronskian positive frequency 
solution of \d 2 + w 2 ] $ p = 0. Evaluated at late time 
(p/a <C Hp), one obtains the standard expression 



' p 



Hp 
2p 3 



(8) 



where H p is the value of H when the p-mode exits the 
Hubble scale. 

The goal of this paper is to compute (numerically 
where necessary) the modifications of this power spec- 
trum due to interactions with some additional field VP 
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inducing dissipative effects as parametrized in eq. {2}. 
To this end, we use the modei introduced in [14;] whose 
essential feature is that + Si nt , the action of \& plus 
that governing the interactions, breaks Lorentz in- 
variance in the UV sector. As a result, similarly to the 
dispersive models of eq. (pj, the propagation of $ re- 
mains unaffected in the low-energy sector, whereas it is 
no longer Lorentz invariant in the UV, even in the vac- 
uum. 

There is of course a lot of freedom to choose S-q, and 
S'int. But as explained in the Introduction, our aim is to 
obtain simple equations for the effective propagation of 
$ after having traced over 4". From this point of view, 

is introduced only to give rise to dissipative effects 
while preserving unitarity. Since Gaussian models are 
the simplest and yet do the job, we work with quadratic 
actions. For further discussion concerning the general 
character of these actions, see [lH and Sec. Ill Fl below. 

Because of Gaussianity, the action still splits into dis- 
connected sectors: 

S=^Jd 3 p (S* (p) + (p) + S int (p)) ■ (9) 
The action 5$ is given in eq. ((H), and we use 

(p) = J d* J dk *t ifc _ (7rA fc) 2 ] * p , fc , (io) 

Si„t (P) = / d V g p ( V ) J dk ($p^*J ifc + h.c.) . (11) 

In the action (fT0")l . the proper frequency of the ^ is 
dimensionalized by A which is the only constant (proper) 
scale. These frequencies are chosen to remain constant as 
the universe expands because this guarantees that the 'J'fc 
are not excited by the cosmological expansion. Moreover, 
the carry no spatial momentum. They are thus at rest 
with respect to the cosmological frame. 

In eq. (fTTj) , the coupling is bilinear, so that the model 
is indeed Gaussian. Hence, we can integrate out the 
^Pfc. Moreover, the index k is chosen to be continuous 
so that the Poincare recurrence time of the system is in- 
finite in order to effectively obtain dissipation 25]. The 
time dependent coupling g p will be chosen so as to pro- 
duce the desired dependence on P/A in the decay rate T, 
see Sec. Ill Bl In order to obtain an effective equation of 
motion of $ which is local in time, see eq. (fTg|) below, we 
have used a derivative coupling. 

B. The effective equation of motion of <3? 

Since our model is Gaussian, all equations of motion 
are linear in the field amplitude. Hence, they can be 
treated as equations for the field operators in the Heisen- 
berg picture. 

The equation for ^>k is 



where D, k = irA \ k\. Its general solution is 
*p,*(*) = * P>fe (*)- /dt' Gt(t,t')d t ,(g p <S> p ) . (13) 

Here, *f? p k obeys the homogeneous equation, and the re- 
tarded Green function is given as 



* _ . dn e - <n ('-'') 

u k V J > t ) - I — 7v2 



2tt Q.% - ft 2 - iefl 



(14) 



Similarly, the equation for $ reads 

[d 2 v + u 2 p (77)] $ P = J dk g p d n y p>k . (15) 
Inserting eq. (fl~3|) into the r.h.s. yields 

[£>l + (»?)] % = Jdk g P d^° Ptk - 9 fd v ( 5p $p) , 



where we have used 



J dkd t ,Gt{t,t') = - 



S(t- t') 
A 



(16) 



(17) 



It is this equation which has motivated our choice of the 
action (fit))) and (fTTj) . Indeed, in general one would have 
obtained a non-local equation, whereas here, the effective 
equation of motion of $ p is simply 

[d 2 n + 2 lp 8 v +u 2 p + d nlp ] % = g p J dk d^ k . (18) 

The term 2j p d v gives rise to dissipative effects. They 
are governed by the decay rate (in conformal time) 



7p0?) 



{9 P {v)Y 
2A 



(19) 



The time dependence of the coupling g p is fixed by the 
following conditions. We first demand that the scale A 
be the proper energy at which interactions between $ 
and appear, irrespective of the comoving momentum 
label p. We also impose g p — > for low momenta P <C A, 
so that the $ p decouple from <J and propagate freely. 
This is required by particle physics observations which 
put severe constraints on possible violations of Lorentz 
invariance, see e.g. fifty . 

In analogy to eq. (UJ, in order to cover the general 
case we classify dissipative effects according to the lowest 
order of P/A: 



7p _ £ 
p P 



l + O 



(20) 



[9; 



(12) 



The first series coefficient can always be set to unity by 
a redefinition of A. From this equation one can already 
conclude that in cosmology, as the proper momentum 
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P = p/a redshifts, the modes go from a strongly dissipa- 
tive regime Y/P — 0(1) for P > A, to an underdamped 
regime where Y/P <C 1. Using eq. (fl"9j) . we see that the 
n-th coupling function should be taken as 



9 P 




n/2 



(21) 



In this we follow the same approach as previously 
employed in dispersive models. First, we replace the 
relativistic relation by an effective equation of motion, 
cq. (|18[) . where the dissipative effects are chosen from 
the outset, and second, we determine the modifications 
of the power spectrum induced by this replacement. In 
this paper we have no ambition to put forward "privi- 
leged" (or "inspired") dispersive/dissipative models that 
could be derived from first principles. Yet another way 
to position this approach is to state that we follow a 
bottom- up rather than a top-down route to new physics. 



C. The power spectrum in dissipative settings 

The general solution of eq. (|18[) is 

*p(»?) = + Jdr/G*(r,,rf) Jdkg p d v ^ k , (22) 

where $p(f?) and G* are the homogeneous solution and 
the retarded Green function, respectively . 
The homogeneous solution <&p decays as 



$£(77) oc exp 



>lin 



(23) 



where the initial time r]i n fixes the moment when $ and 
^ start to interact. Since we do not want to fine-tune 
the number of e-folds, we assume that ryin is located deep 
in the ultraviolet regime 



Pn 

A 



P 



a(?7in)A 



> 1. 



(24) 



In this case, the $p does not contribute to any observable 
at late time, implying that only the state of the is 
relevant [26| . Therefore, the power spectrum of super- 
horizon modes is insensitive to the initial state of <f> p . 

Let us establish this important property in more detail. 
If the initial state factorizes, i.e. \& p k and $p are not 
initially correlated, the anti-commutator of $ reads 

({%(v),%'(v')}) = ({<S> d p (v), *£(»/)}) 
rn' 

dm / dr? 2 G ! *(r7,7?i)G*,(r?',?72)A r p,p'('7i,?72), 

riin Jvin 

(25) 

where we introduced the so-called noise kernel 
AA p ,p, (771,772) = 5 3 (p-p')N p (771,772) = 

g P (m)g P (m) [ dkdk'({d m *l k ,d m V p \ k ,}) , (26) 



whose properties will be specified below. We have simply 
assumed that the state of ^ k is homogeneous. Because of 
the decay given in eq. ([23]) , we immediately conclude that 
the first term in eq. (|25p will be exponentially damped. 
Hence, the anti-commutator, and thus the power spec- 
trum ([6]), are entirely given by the term which is driven 
by the noise kernel. 

In other words, if inflation lasts long enough and if 
dissipation is sufficiently efficient in the UV (this require- 
ment will be discussed in more detail below), the equation 
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P p = lim 

» )fin -*o- 2a z (7j fin ) 



rVlin rVtix, 

\ I d?7i / d?72 G*(?7fin, Vi) Gp(v&n, m) Np{rn,r) 2 ) , 



'/i. 



mi 



(27) 



is exact, and we may furthermore take 77i n to — oo. This 
equation replaces the usual expression of eq. 0, gov- 
erned by the norm of the free mode <I>p a , and valid both 
for relativistic and modified dispersion relations. 

In view of eq. (|2T|) . we see that dissipation affects the 
structure of the equations much more profoundly than 
dispersion does. We also understand that the introduc- 
tion of ^ could not have been avoided, since 'J determines 
both the noise kernel N p (through its anti-commutator) 
and the decay rate 7 P (through its retarded Green func- 
tion, see eqs. (fiUl UM). These must be related to each 
other by a fluctuation-dissipation relation, see [l4[ for a 
brief review in the present context. This explains why, 
unlike dispersion, one cannot treat dissipative effects by 
simply introducing an imaginary term in the dispersion 
relation. 

The remarkable property of dissipation when it is intro- 
duced by coupling to some dynamical degrees of freedom 
is that P p is independent of all their properties if H p -C A 
and if they are in their ground state. Moreover, in this 
case, as we will show, the spectral power (l27l) agrees with 
the standard value given by eq. (J5J . 



D. The noise kernel 

The definition of the model is complete once we specify 
the state of the VPfc. For simplicity, we only consider 
thermal states. Recall that the proper frequencies of ^>k 
are time independent so that the proper Hamiltonian of 
has stationary eigenstates. 

Then, at temperature T, the noise kernel (j26|) is 

N p = g P (771) g p (772) a (771) a (772) 
2T 

x —dt 1 fx/bh(nT{ti-t 2 )) ■ (28) 

This directly follows from the fact that the (free) fields 
can be decomposed as 

*p,*(*) = 6p,fctffc(*)+&! Pj - fc (29) 
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where a Pi fc satisfy canonical commutation relations and 
where the isotropic mode functions 



ipk = 



1 



^—iQkt 



(30) 



have unit Wronskian WfV'fc] = 2lm(ipkdtipl) = 1. We 
have also used the fact that in the thermal states one has 



(31) 



At high temperature, as usual, the kernel becomes local 
lim N p = 4T lp ( m ) a(rn) Sfa - m) > (32) 

T — >oo 



where we used 8{t\ — t 2 ) = S(r)i — r] 2 )/a(rji) and eq. (|T9 
At zero temperature (vacuum) instead, one gets 



4 / PV 

lim i N p = -Jj P {mhp(m)a(m)a(m)dt l7 -. — -t, 

1 ->0 7T V [ti — 12) 

where the singular behavior should be interpreted as the 
derivative of the Cauchy principal value. 



when the constant Wronskian of \p was chosen to unity. 
It should be noticed that only the decay accumulated 
from 77 to 7/ appears in Gf. The fact that r/o must drop 
out can be seen from eq. (|3"4"|) . In fact, the second equal- 
ity replaces the equal time commutation relation in the 
presence of interactions, see [l4[ for further details. 

Returning to eq. (12"?| , the presence of the two functions 
X p evaluated both until ?7fi n = limits the past history 
that is relevant to the power spectrum of super-horizon 
modes. To characterize this relevant domain, we define 
the time rj* by the moment where 



1 P (r? fin = 0, 77*) = 1 . 



(40) 



Times earlier as rj* play no significant role in the power 
spectrum. In other words, T p can be considered as an 
"optical depth" . 

Having established these features, we can now explain 
why, if H <C A, any dissipative model exhibiting dissi- 
pation above A behaves as if it belonged to the class of 
models we just considered. 



E. Retarded Green function 



To compute the power spectrum (|27|). we need the re- 
tarded Green function of eq. (fT8|) . It satisfies the bound- 
ary conditions 

G£(t/ = tj)=0, d v G%,= v = l. (34) 

Therefore, it can be written as 



-..|.,„ _ u,.. „n 21m (p„ (7/) (r/0) 



where the mode function tp p (rf) may be any homogeneous 
solution of eq. (fT5)) that has a nondegenerate Wronskian 
W[<p p }. 

We introduce the function 



2p iv,m) 



1p iv') dr?' , 



(36) 



which gives the amount of dissipation from 770 to 77. It 
will play a crucial role in what follows. Using it, we can 
get rid of the friction term in eq. I|18p , by writing 

<P P (V) = e- T ^Xp(v) ■ (37) 
Indeed, \p obeys 

[d^+^ p (v)-ll(v)]x P (v) = 0. (38) 

Taking into account the time dependence of the Wron- 
skian of ip p , eq. (|35|) can be rewritten as 



G* (77,77') = -20(77 - rf)Tm [x P (v)x* P (v')] e 



(39) 



F. General properties of dissipative models 

The models we studied are based on several simplifying 
assumptions. First, they are Gaussian; second, the fre- 
quency of the is constant; and third, a derivative was 
introduced in the action Si nt in order to get a local equa- 
tion for <!>. Nevertheless, the features we obtained are 
more general: they will be found in all dissipative mod- 
els respecting minimal assumptions that we now clarify. 

Before listing these conditions, it should be noticed 
that when dealing with nonlinear interactions, it is no 
longer convenient to work with the mode operator $ p as 
we just did. Instead, it is appropriate to study the effec- 
tive evolution in terms of the two-point correlation func- 
tions of <i>p, see |l3 |. In particular, it can be shown that 
the expectation value of the anti-commutator of the 
l.h.s. of eq. (|25|) . always obeys a linear integro-differential 
equation with a source [171 ]. This also applies to non- 
derivative, bilinear couplings, so that the following dis- 
cussion includes both cases. 

Adopting this language, we can transpose the two con- 
ditions we used in Section III CI First, the dissipative 
effects should be strong enough so as to erase the con- 
tribution of the homogeneous solution of this integro- 
differential equation. In this we recover, in the lan- 
guage of two-point correlation functions, the neglect of 
the homogeneous solution $ p of eq. (fig)) . Second, the 
state of the entire system must be spatially homoge- 
neous. When both conditions are met, the expectation 
value of the anti-commutator of the Fourier mode $ p is 
driven by a p-dependent (c-number) source through the 
above-mentioned integro-differential equation. This im- 
plies that the power spectrum will be given by eq. (|27|) 
in any (unitary) dissipative model, Gaussian or not. In 
other words, the power spectrum is always governed by 
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a p-dependent kernel N p and a retarded Green function 
p ' 

Let us begin with the kernel N p , which encodes the 
properties of the state of the system. In Gaussian models, 
it is simply given by the expectation value of the anti- 
commutator of the r.h.s. of eq. (fT8|) . In non-Gaussian 
models, it must be computed order by order in a loop 
expansion. This calculation might turn out to be diffi- 
cult, but (in renormalizable theories) N p is a well defined 
kernel which is given by the real part of the (renormal- 
ized and time-ordered) self-energy of $ p (see for instance 
Appendix B in [HI). Therefore, when N p has been com- 
puted, it will "drive" the power spectrum as indicated in 
eq. ([27]). 

Let us briefly discuss the modifications one encoun- 
ters when the proper frequencies f2fc(£) of the de- 
pend on time. In this case, their state will be para- 
metrically excited. However, if the variation of 0/j(£) 
is slow enough, this amplification will be exponentially 
suppressed by virtue of the adiabatic theorem. (A simi- 
lar situation is expected when dealing with non-Gaussian 
models.) Then, if the are initially in (or close to) their 
ground state, N p will essentially be the noise of the adi- 
abatic vacuum of the \I/fc . 

Let us now turn to the effects of dissipation. In a gen- 
eral model, one would lose the local character of eq. (jTHJ). 
However, in all models (Gaussian or not), the retarded 
Green function of $ obeys a linear integro-diffcrcntial 
equation of the form 

[dl + LJ 2 p ] G*(ti, rf) + J V d m V p (rj, m ) Gp (771 , rf) 

= 6(r)-ri'), (41) 

where the non-local kernel T> p generalizes what we had 
in eq. (fT5)) in that, when V p — d m S(r) — r/i)2j p , one 
recovers the usual odd term of that equation. The kernel 
T> p is antisymmetric in the exchange of its arguments 
and describes dissipation. It is related to the imaginary 
part of the time-ordered self-energy, as N p was related 
to the real part, and is therefore also well-defined and 
computable, at least perturbatively. Moreover, as for 
N p , if the state of the evolves adiabatically, T> p is the 
dissipation kernel in the adiabatic vacuum. 

Concerning the power spectrum, we saw in Sec. Ill El 
that in expanding universes with H/A <C 1, only the 
evolution in the underdamped regime is relevant. In 
this low-energy, weakly dissipative regime, the non-local 
equation (|41[) can be approximated by a local one (i.e. 
similar to eq. (fl"gj) with an effective damping rate j c g) 
provided the characteristic comoving time describing the 
retardation effects of T> p is much smaller than oj" 1 , or 
equivalently that the corresponding cosmological time is 
much smaller than ff~ 1 .This approximation is similar to 
the diffusion approximation in kinetic theory. In this 
case, we can approximate 

/ dqi V p {t), 7?i) <p(rn) ~ 2-f cS (r})d n (p(r)) , (42) 



where 7 e ff(??) depends in general on 77 and on the state 
of the system. (In the case one would consider the same 
model in Minkowski spacetime, and in its ground state, 
the above equation is easily obtained in the frequency 
representation by performing a Taylor expansion in the 
frequency and truncating at first order. In that case, 7 g 
would be constant. In an expanding universe, it becomes 
time dependent through the scale factor a(t).) 

In conclusion, a general <E»- 1 3/ model where (i) H/A <C 1, 
(ii) the state of ^5>k evolves adiabatically, (iii) the charac- 
teristic time of V p is much smaller than H~ x so that 
eq. (1421) holds, will give the same power spectrum as 
that of the corresponding simplified model governed by 
eqs. (jMTTj) with the coupling g matching the effective 
decay rate 7 e ff through eq. (TTJ 



III. ANALYTICAL TREATMENT 

We present some analytical expressions for the power 
spectrum which will facilitate the interpretation of the 
numerical results. They are valid for H <C A and in the 
slow-roll regime. 



A. More properties of the retarded Green function 

As we are interested in the power spectrum of super- 
horizon modes, we can make a first approximation by 
factorizing the growing mode of Xp(v) : 



iH p a 



(43) 



By eq. (|20|) . this is a solution of eq. (|38|) for |pr?| <C 1. 
Hence, for rjan — > 0~ we have 



(-»?) 



2H p a(77 fln ) 



x Re[x P (f7)] e- 1 *^^ . (44) 



G p (as function of rf) oscillates with a slowly varying 
envelope. 

We now give an analytic approximation to eq. (|44p 
valid in the slow-roll approximation and in the case of 
scale separation: 



H 

e,-«l. 



(45) 



Let us first give an approximation for the envelope, given 
by the exponential in eq. (|4"4"|) . We take 7 P to be of the 
form ([20]) . i.e. j p oc a~ n for modes below the UV scale A. 
The term coming from the upper bound is then negligible 
since J p (t)) — * for 77—5-0. Hence during inflation, the 
integral is dominated by the lower bound, and we may 
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thus estimate 



a(i)fin) 

1 



Jp_ 

Ha 2 
1p(v) 



da 



n + 1 i? (77)0(77) 



(46) 



where we have used the slow-roll approximation. 

Let us now consider the term Re [x P (77)] ■ The equa- 
tion of motion, eq. (|38p . has an oscillating solution in- 



side the horizon. However, the term 7^ in the effective 
squared frequency, i.e. the frequency shift due to dissipa- 
tion, introduces some non-adiabaticity to the evolution of 
the mode close to the time when it leaves the ultraviolet 
regime. This implies that the retarded Green function re- 
ceives non-adiabatic corrections for very early times and 
may even stop oscillating in the case where there is an 
overdamped regime (7^ > w 2 ) in the UV. However, in 
the case of scale separation this ultraviolet behavior of 
the mode function occurs only where the envelope is ex- 
ponentially small by a factor e~°^ A ^ H \ In other words, 
the dispersive effects induced by dissipation are damped 
by dissipation itself (see also fig. [5]) . 

Finally, notice that the noise kernel ([28]) is propor- 
tional to 7 P and therefore vanishes as 771,772 — > — . In 
conclusion, the double integral (j2"T|) takes its value in the 
vicinity of 771, 772 ~ 77* defined by eq. ([40)) . If H -C A, one 
finds that XpiVp) is wei l inside the horizon and at the 
same time sufficiently below the UV regime, such that it 
is justified to use a unit Wronskian free oscillator in place 
of Xp m order to estimate the power spectrum. In other 
words, as an approximation, we set 



Re [xp{v ~ Vp)] 



1 



cos(p?7) , 



(47) 



and thus 



Gp(?7fin^0 ,77^77* 



x cos(pr])e "+ 1 HiWi) (48) 



when we evaluate the double integral. It should be 
stressed that this is only done to get an analytic estimate, 
and that no approximation is required for the numerical 
treatment since we can always solve eq. (|38|) numerically. 



B. General properties of the power spectrum 

Inserting eq. (I48[) into eq. ([2"?)) we have 

,0 ,0 1 

' P p = ' P p / dr li / d V2 -cos(p77i)cos(p77 2 ) 
J —00 J— 00 P 

1 / -lp(vi) 1 Jp(V2) \ 



x e 



The factor a 2 (f]fi n ) coming from the rescaling of cf> is 
compensated by the two factors of a(?7fi n ) generated by 



the growing modes of Xp- The integrand is now indepen- 
dent of 77fj n . 

In the vacuum, i.e. at vanishing t emp erature, it has 
been shown by general arguments [14l . Il9j ] that the power 
spectrum agrees with the standard prediction JHJ for 
A ^> H. In other words, the double integral on the 
r.h.s. of eq. (|49|) evaluates to unity in this double limit. 
This can be shown by a lengthy calculation which will be 
omitted here. We only state that a naive analytic esti- 
mate can be found [2(| for the magnitude of the leading 
order modification: 



with SV V 



dlnSVp 
d\n{H/K) 



1 p p • 



O 



A" 



(50) 



C. High temperature limit 

In keeping with the general character of our model, we 
do not specify a physical motivation for the case of finite 
temperature, where $ acts as a "heat bath". However, 
it serves to illustrate the mechanism how the state of the 
adiabatic modes is dynamically determined by the state 
of the In particular, if the latter are in a thermal 
state, $ thermalizes through the interaction. Further- 
more, we establish that the thermal excitations of <& p 
are effectively populated at 77 — rj* defined by eq. PP]) . 

Let us consider the high temperature limit of eq. (|49D , 
i.e. we insert the limiting expression ([5^]) for the noise 
kernel and set the (5-function against one of the integra- 
tions. We then have 



V P (T^^) = V° p 

ip(m) 



x 4T 



drj] 



2 ~fp(vi) 



a(?7i)cos (prji)e stt»»Cii) . (51) 



Assuming that A 3> H, we may replace the squared co- 
sine by its average value 1/2. Since j p oc a~ n in the rele- 
vant domain, the remaining integral is now essentially a 
representation of the Gamma function. It evaluates to 



Pp(T 



°°) = V p x t(t?- 
p A\Hn 



1 



1 



(52) 

If we compare this result to the high temperature limit 
of a thermal power spectrum 



V p (T) = P°x coth 



2T 



(53) 



where Q* denotes the proper frequency at the instant 
when the occupation numbers are fixed, we find that 



n* 
T 



(n+i) 



1 H 



A 



(54) 



8 



One verifies that f2* so denned coincides with the proper 
frequency at 77*: 

n* * ^ . (55) 

Note that it is below the scale A (see eq. (|54[l ). in agree- 
ment with the results of Sec. IIII Al 

In conclusion, the state of <j> is inherited from that of 
iff at the time 77*. 



6. The power spectrum is calculated from eq. (|56|) . 
The derivatives and the double integral are evalu- 
ated numerically. 

To proceed we first need to address two technical is- 
sues. The first concerns the singular behavior of the noise 
kernel, the second the asymptotic behavior of 7 P in the 
ultraviolet. 



IV. NUMERICAL ANALYSIS 

We present a numerical scheme that solves the dou- 
ble time integral of eq. (j2"T|) by means of Monte Carlo 
integration. 

A. The numerical scheme 

The procedure includes the following steps: 

1. We impose the inflationary background. Specifi- 
cally, de Sitter space and power law inflation are 
considered. 

2. We specify the "relative" decay rate j p /p in terms 
of a function of P/A, see eq. (|2H)) . 

3. The equation of motion eq. for \p (v) is solved 
numerically. 

4. The integral X v (77, 77') of eq. ([36]) governing the 
amount of dissipation is computed numerically. 

5. The retarded Green function is constructed from 
these numerical solutions, cf. eq. (|39[) . 



B. Handling the singular behavior of N 

The noise kernel (cf. eq. ([28]) ) is singular for equal times 
and should be interpreted as a Cauchy principal value. 
Unfortunately, numerical integrators are generically inca- 
pable of calculating principal values; they usually fail to 
achieve convergence within a finite number of integrand 
evaluations. However in our case, the integrand can be 
rewritten in a regular form by a convenient change of 
variables. 

We first note that the integrand of the double integral 
in eq. (|27p is symmetric under the exchange rji «-> 7/2. 
Thus, if we make a change of variables to £ = 7/1 + 7/2, 
£ = 7/1 — 772, the integrand will be symmetric around 
£ = 0, where it has a double pole. Following the general 
techniques in the calculus of generalized functions [2l| . 
the integrand may be regularized by multiplication with 
— £ 9j£ = 1 and performing an integration by parts in 
£. The derivative now acts on a regular (the £ 2 cancels 
the double pole) and symmetric function. It thus gives 
an odd function in £ such that the remaining singular 
factor in £ _1 may be lifted. This way we obtain a regular 
integrand that is well-behaved within the entire domain 
of integration. In brief, we have 



/ydCd£ 



2a 2 ( VfiD )JJ 2 e * 



G v 77fi n , 



N r . 



2 ' 2 



r 1 - 



(56) 



C. The UV behavior of the coupling function 

For numerical integration, the range of 77 has to be 
truncated somewhere in the remote past. In order to 
guarantee a safe truncation, we impose that the inte- 
grand drops off exponentially. Then the cutoff can be 
chosen in such a way that the truncation error is negligi- 
ble w.r.t. the numerical value of the integral. 

The exponential behavior of the integrand is achieved 
by a suitable choice of the decay rate j p . One might 



be tempted to say that any positive 7 P gives rise to an 
exponential behavior, simply given by T p . However, the 
effective frequency of the mode functions depends on the 
damping rate as well, cf. eq. (f3"5)) . If j p is not bounded 
from above, we may run into some pathologies due to an 
unbounded ovcrdamping in the UV. The reason is that 
dissipation is less effective in an overdamped situation. 
To see this, let us consider a classical oscillator with con- 
stant frequency lu and damping rate 7 in the overdamped 
regime, i.e. 7 > u>. It has two decaying modes with the 
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decay rates 7 ± y/j 2 — u> 2 . When j 2 3> lo 2 , the slowly 
decaying mode is thus cx exp[— ^- (77 — 770)]. 

Returning to our model, assuming j p /p — (p/aA) n 
without higher order terms, the growing WKB solution 
of eq. (|38|) for j p ^> oj p would lead to an overall asymp- 

totic behavior cx exp(— H 27 Fn 7 ] ^')- ^ s tne m tegrand 
drops off like a n in the remote past, the integral is gen- 
erally finite in the limit rj ln —> — 00. In this case nothing 
guarantees that one can disregard the damped initial cor- 
relator in eq. (|2"5")1 and that the power spectrum does not 
depend on the initial state of $ at r/j n . 

However, when A 3> H, the residual contribution 
of the decaying term in eq. (|25|) is strongly suppressed 
(0 [exp (— A/i?)]) because of the dissipation between A- 
crossing and horizon exit, where the mode is in the un- 
derdamped regime. But if we want to neglect the con- 
tribution of this decaying term for all values of the ratio 
A/H, we cannot work with j p /p = (P/A)" (nor with a 
polynomial of finite order) since it is not bounded from 
above. We will instead work with a decay rate that sat- 
urates in the UV. We choose 



If 
P 



— = K tanh 



((- 

VVaA 



(57) 



as a simple realization of this property. The new param- 
eter k was introduced such that it only appears in the 
subleading terms and that the decay rate saturates in 
the UV when it reaches the value j p /p — k. 




in units of H] 



FIG. 1: Power spectrum V p (in units of the standard power 
V°) for n = 2 and T = in de Sitter space, as a function 
of A/H and for various choices of k. For A 3> H, the power 
spectrum asymptotes from below to the standard value ir- 
respectively of the value of k. For A <C H, the spectrum 
becomes flat again but it is lowered by a K-dependent factor. 
The region A < H is non universal. 



D. Numerical results 

1. Dependence on H/A in the vacuum (T — 0) 

Let us now consider the zero temperature limit, i.e. ^ 
is in its vacuum state. On physical grounds, we expect 
that this will be the relevant case if a fundamental theory 
gives rise to dissipative effects in the UV sector during 
inflation. 

The modified power spectrum is computed for A rang- 
ing from A 3> H down to A <§C H and is compared to the 
standard prediction f8j, both in de Sitter space (figs. [T] 
and [2]), and in power law inflation (fig. [21) ■ 

Let us first discuss the power spectrum in de Sitter 
space as a function of A, fig. [TJ coming from the high 
values. Note that in de Sitter, A/H is time independent, 
so the power spectrum is scale invariant and the value 
of A only affects the normalization of the spectrum. As 
expected from the analytical results, for A 3> H the stan- 
dard power is recovered, independently of k. In this we 
corroborate the robustness of the power spectrum when 
the initial state is the adiabatic vacuum [12j (see also 
below in this section). As A approaches H, the power 
spectrum is modified in a non-universal way which de- 
pends on all model parameters. The flattening of the 
curves for A <C H has to be attributed to the fact that 
our particular choice of j p (cf. eq. (|57|l ) saturates to a 



constant in the UV and thus becomes independent of A 
inside the horizon. 

Having established the robustness in the regime 
H/A <C 1, the signature of dissipation is contained in 
the first deviation with respect to the standard result. 
As anticipated in eq. ([SH]) . in the limit H/A — > 0, the 
deviation behaves as 



V* n { A 



(58) 



where the constant S n depends only on n. This has been 
verified for values of 1 < n < 2.5, as defined in eq. (|57| . 
and we conjecture that it is valid for any power. For all 
these values, 8 n is negative which means that the spectral 
power is reduced w.r.t. the standard value. We were not 
able to probe higher values of n because of the numerical 
difficulties to follow the sharp decrease of the modifica- 
tions. However, preliminary results indicate that (as for 
dispersive models [13j) 5 n changes sign for n — 3, which 
means that for n > 3, dissipation leads to an increase of 
the power. At present we have no explanation for this 
unexpected result. 

In fig. the difference 8V P = V v — V p is plotted for 
n = 2 and various values of k using a logarithmic scale. 
It shows that the power spectrum becomes insensitive to 
the next-to-leading order terms in the expansion of j p , 
see eq. (|2T)j) . 
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Comoving momentum p [arbitrary units] 



FIG. 2: Difference between the modified (n — 2, T = 0) and 
the standard power spectrum in de Sitter space for A > H. 
The straight line indicates the first order deviation (which is 
independent of k) with a slope equal to n. 



Let us now turn to power law inflation figs. [3] and [4] 
Fig. [5] indicates how the slow-roll parameter e impinges 
on the features of the modification SV p . Surprisingly we 
find no dependence of S n on the slow-roll parameter, or if 
present, it must be extremely mild. This implies that the 
power spectrum from power law inflation, as a function 
of p (fig. [3]), is the power spectrum in de Sitter, as a 
function of Hp/A, where H p is the value of H when the 
p mode exits the horizon, up to 0(j 2 /p 2 ). 

Note that the modification of the power spectrum due 
to dissipative effects can be reinterpreted in terms of a 
(scale dependent) spectral index, which is given by the 
slope in fig. [21 To lowest order in e and H p /A, using eq. 
(1551). we find 



dln(p 3 V p ) 
alhip 



2e 



(59) 



The first term accounts for the slow-roll evolution of the 
standard power spectrum, whereas the second term is 
due to dissipative effects. It is proportional to e and 
suppressed by (H p /A) n <C 1. 

For n = 2, the negative sign of 5 n gives rise to a down- 
ward modification of the spectral power. Moreover, the 
effect is enhanced at large scales due to the p-dependence 
of H p / A. The effect tends to increase the spectral index 
(i.e. increase the slope at any given point of the curve), 
and the running of the spectral index should indicate a 
concave spectrum (the curve is bent downwards towards 
large scales). 



FIG. 3: Power spectra for power law inflation (e = 0.2). n — 
2 and T — 0. The continuous line indicates the standard 
power p 3 Vp. The vertical line corresponds to the comoving 
momentum p where H p = A. For H p <C A the power spectrum 
asymptotes to the standard value, whereas in the other limit 
it is again sensitive to the UV. 



2. Dissipation and non-adiabaticity 

We now discuss the interplay between dissipation and 
non-adiabaticity in more detail. Following what was done 
with dispersive models [l3| , we plot the function 



a 



d v LU e f{ 



as a function of x = —pn in fig. O where 



Weft 



1p ' 



(60) 



(61) 



is the effective (conformal) frequency of the modes x P , 
see eq. ([55)1 . We see two regions of non-adiabaticity, the 
usual one at the horizon x < 1, and a new UV feature, 
inside the horizon, which is due to the dispersive effects 
induced by the dissipation and governed by -f 2 . We also 
plot the exponential factor exp[— 1(0, rj)], see eq. (|36l) . 
appearing in the retarded Green function. We see that 
for sufficiently large A, the UV feature in o~(x) lies in the 
tail of the exponential. The non-adiabatic effects caused 
by dispersion are then completely masked by dissipation. 

As A decreases, the UV feature in a{x) begins to over- 
lap with the exponential. For A ~ H, it is well inside 
the region where the exponential is still 0(1). In that 
case one expects significant modifications of the power 
spectrum. 
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A [in units of H] 



FIG. 4: Difference between the modified (n = 2, n = 2, T — 
0) and the standard power spectrum for power law inflation 
and with A > H p . Only the subleading modifications are 
sensitive to the value of the slow-roll parameter e. 

3. Dependence on T/A 

In this section we consider a de Sitter background and 
set n — 2. In fig. [6] we plot the power spectrum as a func- 
tion of T in a case H <C A. For T < H it asymptotes 
to the standard power T^, while for T > A it converges 
to the analytic result ([52]) . We verify that this curve is 
fitted at better than 1% by the power spectrum obtained 
in non-dissipative and relativistic settings when the "ini- 
tial" state of <E> is enforced to be that thermal state at 
r] = rj* see eq. (l53l) . 



V. CONCLUSION 

Let us sum up our main results. First, in spite of the 
strong dissipative effects encountered in the early mode 
propagation, the predictions for the power spectrum con- 
verge to the standard ones in the case of scale separation, 
Hp -C A, and if the environment field \& is in the ground 
state. The power spectrum is hence a robust observ- 
able with respect to high-energy dissipative effects under 
these conditions. In this regard, dissipative models do 
not differ from dispersive ones [12j. 

Second, the leading deviation of the power spectrum 
induced by dissipation is linear in the relative decay rate 
1p/p evaluated at horizon crossing, see eq. (|58p . The sig- 
natures of dissipation therefore do not oscillate. In this 
regard as well, dissipative models behave like dispersive 




x = —prj 



FIG. 5: The non-adiabaticity coefficient a and the exponen- 
tial ([361) for k = 2 and k = 0.5. Upper panel: A/H = 10 3 , 
i.e. the scales are clearly separated. Lower panel: A/H — 10, 
i.e. weak scale separation. The background is de Sitter, and 
n — 2. [a was enhanced by a factor of A/H for a clearer 
presentation.) 



ones [13( • In the region of parameter space we succeeded 
to explore, i.e. the power n < 2.5, the deviations are neg- 
ative, giving rise to a suppression of the power spectrum 
with respect to the standard one. 

Third, we have verified that the deviations in power 
law inflation essentially behave as those evaluated in de 
Sitter space, with H replaced by H p evaluated at hori- 
zon exit. As a direct consequence, any observable effects 
of high-energy dissipation will be more pronounced at 
the largest accessible scales, corresponding to the largest 
Hp/A. 

Fourth, in the case of H p <C A, we showed how dissipa- 
tion, via the function X p , sets the time when the "initial" 
state of the $ p mode is effectively set. It is given by 
eq. (|4"0|) where P* = p/a(r)*) is at an intermediate scale 
between H p and A. In addition, the state of <j> p coincides 
with that of the degrees of freedom causing the dissipa- 
tion. Thus, the properties of the state of ^>k pass on to 
that of <I>p at 77* . 

Fifth, in Sec. Ill F[ we gave (sufficient) conditions under 
which a model exhibiting dissipation above the scale A 
can be well approximated by a Gaussian model, in the 
sense that both models predict the same power spectrum. 

Our analysis could be extended in two directions. 
First, we considered a Gaussian model. If we relax this 
hypothesis, since dissipation grows with the coupling, it 
will be interesting to investigate the combined effects of 
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Environment temperature T [in units of H] 



FIG. 6: Power spectrum as a function of the temperature T, 
computed for a de Sitter background with A = 100 H and n = 
2. The result is therefore independent of k. For T <C H , V v 
asymptotes to Pp (horizontal dashed line), whereas for T ^> 
H it asymptotes to the high temperature limit, eq. (|52[) (the 
oblique dashed line). The numerical result neatly reproduces 
the power spectrum of a thermal state (cf. eq. (|53f> and below) 
whose temperature is specified at rj*. 



dissipation and non-Gaussianities. Second, we calculated 
the power spectrum of a test field propagating on an 
inflationary background. It is a challenge to construct a 
realistic model of inflation displaying dissipation in the 
UV sector. 

We can nevertheless make the following observations. 
At the linearized level, the spectrum of scalar metric per- 
turbations £ is related to that of our scalar field </> by 
V v = ^T V i & This im plies that the relative modifi- 
cation SV^ [Vi due to some dissipative effects is the same 
as the one of our scalar field. Therefore, if \& couples iden- 
tically to scalar and tensor perturbations, the S/T ratio 
should not be changed at first order. Then the absence 
of features in the power spectra in the regime H/A -C 1 
would prevent us from disentangling the new physics 
from a simple shift of the inflaton potential, adding yet 
another ambiguity to the program of the inflaton poten- 
tial reconstruction. 
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